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ABSTRACT 


bi 


A computer-aided flow visualization process has 
been developed to analyze video images acquired 
from rotating and translating light sheet visual- 
ization systems. The computer process integrates 
a mathematical model for image reconstruction, 
advanced computer graphics concepts, and digi- 
tal image processing to provide a quantitative and 
visual analysis capability. The image reconstruc- 
tion model, based on photogrammetry, uses knowl- 
edge of the camera and light sheet locations and 
orientations to project two-dimensional light sheet 
video images into three-dimensional space. A so- 
phisticated computer visualization package, com- 
monly used to analyze computational fluid dynam- 
ics (CFD) data sets, was chosen to interactively 
display the reconstructed light sheet images, along 
with the numerical surface geometry for the model 
or aircraft under study. A description is provided 
of the photogrammetric reconstruction technique, 
and the image processing and computer graphics 
techniques and equipment. Results of the computer 
aided process applied to both a wind tunnel trans- 
lating light sheet experiment and an in-flight ro- 
tating light sheet experiment are presented. The 
capability to compare reconstructed experimental 
light sheet images and CFD solutions in the same 
graphics environment is also demonstrated. 
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NOMENCLATURE 

A,v,z 

2D 

two-dimensional 

a 

Ax , Ay 

3 D 

three-dimensional 

a; 

affinity distortion terms 

U>, K, <f> 
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coefficients of the general equation of 
the light sheet plane 

camera constant corresponding to the 
lens effective focal length 

input 2D digital images 

output 2D digital images 

radial distortion terms 

pixel gray level value 

3 by 3 rotation matrix 

horizontal scan line position of a pixel 

horizontal scan line position of center of 
brightness 

column position of a pixel 

column position of center of brightness 

lens decentering distortion terms 

parameter in parametric equations of a line 
in 3 D space 

image space coordinates 

photogrammetric principal point 

object space coordinates 

angle of attack, deg 

corrections for removing biases in the 
imaging system 

angles for camera orientation, deg 

light sheet position angle, measured 
counterclockwisg/from fuselage nose, deg 
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Subscripts: 

C camera lens center 

I image point 

0 object point 

S light sheet origin 
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INTRODUCTION 


Flow visualization techniques are important 
tools for the aerodynamicist. The ability to “see” a 
flow can aid in the development of heuristic models 
to describe it. One common technique used at 
Langley Research Center is light sheet illumination 
of particle-seeded flow. A single, stationary light 
sheet allows characteristic structures of the flow to 
be visualized for one slice of the three dimensional 
flow Held. The recording camera is often restricted 
to a fixed location outside the test section for wind 
tunnel applications, or on-board the aircraft for 
in-flight experiments, with an oblique orientation 
to the light sheet. Consequently, the recorded 
images suffer from perspective distortion which 
must be accounted for when making geometric 
measurements. Digital image processing techniques 
are available to extract quantitative data from light 
sheet images, as well as correct for perspective 
distortions. These techniques rely on photographing 
a physical grid of known geometry in the plane of 
the light sheet with the recording camera. While 
this method is relatively simple for experiments 
involving stationary light sheets, it is not practical 
for experiments where the light sheet moves. 

Advances in the design of light sheet systems 
have resulted in greater flexibility of the light sheet 
technique. For example, capabilities now exist to 
control the size and position of light sheets in real- 
time, translate and rotate single sheets through a 
volume, and to synchronize the video with pertinent 
flow oscillations 1 . These systems offer aerodynami- 
cists the opportunity to see multiple slices of the 3 D 
flow field about a model. However, the task of an- 
alyzing and interpreting the recorded 2D images is 
not an easy one. An obvious difficulty lies in trying 
to form a mental image of the 3 D flow field based 
on viewing a sequence of 2D slices through it. Not 
only is visualizing the slices as a volume of data im- 
portant, but also making accurate geometric mea- 
surements of features in the flow field. Clearly, it is 
not feasible to position a physical grid at every light 
sheet plane, especially for experiments where the 
light sheet continuously sweeps or rotates. 

The purpose of this paper is to describe a 
computer-aided process that was developed to en- 
hance the usefulness of 3 D light sheet flow visu- 
alization techniques. The main objective is to re- 
construct the 3 D geometric shape of the visible re- 
gion illuminated by translating and rotating light 
sheets. The computer process is a combination of 
digital image processing and computer graphics 
procedures, and the 3 D reconstruction algorithm 
based on photogrammetry. The process has been 


applied to light sheet images generated at Lang- 
ley Research Center and provides researchers with 
an enhanced analysis capability for understanding 
flow behavior. This paper describes the sequence of 
procedures that have proven to be effective at pro- 
viding quantitative position information for features 
of interest in the flow field, as well as an improved 
qualitative analysis capability. The description of 
the photogrammetric reconstruction technique, the 
image processing and computer graphics techniques 
and equipment, along with results of the computer- 
aided process applied to both a wind tunnel and 
an in-flight flow visualization experiment are pre- 
sented. 


MODEL FOR IMAGE RECONSTRUCTION 


The image reconstruction technique developed 
in this section is based on a simple form of pho- 
togrammetry. It replaces one of the cameras of a 
conventional stereo photogrammetry system with a 
light sheet. The reconstruction model consists of an 
object / image coordinate system transformation, 
the collinearity condition of photogrammetry, and 
the equation of the light sheet plane. A simple set 
of equations relating ideal image and object coordi- 
nates is developed in this section. 

Objects and their images are located with 
respect to an object ( X , Y, Z ) and image (x, y, z) 
coordinate system, respectively. The two systems 
are related by a simple coordinate transformation. 
The origin and orientation of the object system 
is arbitrary. The origin of the image system is 
located at (Xq, Yq, Zq<) in object coordinates. The 
orientation of the image system with respect to the 
object system is given by the rotation angles u>, <f>, k 
via the 3 by 3 rotation matrix M = (mi j) defined 
as 


M = 


cos(k) sin(n) 0 

-sm(s) cos(k) 0 

0 0 1 
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0 10 - 
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The rotation angles are defined as positive for 
counter-clockwise rotation when observed from 
the positive end of an axis looking toward the 
origin. The coordinate transformation from image 
to object system is therefore given by 
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A typical reconstruction geometry is illustrated 
in figure 1. Note that the camera lens is located 
at the origin of the image coordinate system, the 
camera points along the negative z— axis, and the 
x— and y — axis exit to the right and top of the 
camera lens, respectively. The object is located at 
( Xq , Yq, Zq) in object coordinates and its image is 
located at ( xj , yj, c) in image coordinates, where c 
is the effective focal length of the lens which locates 
the image plane. The collinearity condition of 
photogrammetry requires that an object, its image, 
and the lens center lie on a line, as shown in the 
figure. A set of parametric equations for this line 
can be written as 

X = Xc + (Xj- X c)t 

Y = Y c + (Yj- Y c )t (3) 

Z = Z c + ( Zj - z c )t 

where (A/, Yp, Zp) are the object coordinates of the 
image point (xp, yp, c), given by equation 2. These 
equations can represent any point (X, Y, Z ) on the 
line. Since the object point also lies in the plane 
of the light sheet, its coordinates are given by the 
intersection of the line and an equation representing 
the light sheet plane. This intersection process is 
called triangulation. The origin of the light sheet 
is located at (X s , Y s , Z s ), as shown in figure 1. 

In general, the light sheet plane equation can be 
written in the form 

biX + b 2 Y + b 3 Z + 64 = 0. (4) 

By substituting equation 3 into equation 4, t is 
determined to be 

. _ b\Xc + b 2 Yg + b 3 Zg + 64 /gN 

bi(X c - Xj) + b2(Y c - Yj ) + b 3 (Z c - Zj) 

By substituting this value for t back into equation 
3, the unique object coordinates (Xo, Yo, Zq) can 
be evaluated. 

DETERMINING MODEL INPUT PARAMETERS 

In the model described above, the parameters 
for the location and orientation of the camera and 
the light sheet are assumed known. However, these 
parameters need to be determined experimentally. 

In addition, cameras (especially short focal length 
cameras) can exhibit considerable image plane 
distortion which must be accounted for. 

The equations of photogrammetry which model 
the projection of object coordinates (X, Y, Z ) to 
refined image plane coordinates (x, y) of a camera 
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are 


x = x' + Ax = Xp — c 

x mu(Xc ~ X) + mi2(Yc — Y) + mi 3 (Z c - Z) 
_m 3 i(X c - X) + m 3 2 (Y c — Y) + m 33 (Z c - Z) _ 

y = y' + Ay = y p - c 

x m2i(X c - X) + m,22(Yc - Y) + ni2 3 (Z c — Z) 
fn3l(Xc ~ X) + m 3 2(Y c — Y) + m 33 (Z c - Z) 

(6) 

The point (x' , y r ) is the measured image plane lo- 
cation of an object point, (x p , y p ) is the photogram- 
metric principle point, and the A* and Ay terms 
represent functions for modeling the distortions of 
the imaging system. A commonly used set of func- 
tions for A* and Ay is 

Ax = x(k\r 2 + & 2 C 4 + k 3 r 6 ) + pi(r 2 + 2x 2 ) 

+ 2p 2 xy 

Ay = y(kir 2 + k 2 r 4 + k 3 r 6 ) + 2p\xy 
+ P 2 {r 2 + 2 y 2 ) + a\x + a 2 y 

where r 2 = (x' — x p ) 2 + (y' — y p ) 2 = x 2 + y 2 , and 
hi, pi, and a 4 - represent radial, lens decentering, and 
affinity distortion terms, respectively. 

The parameters x p , y p , c, k 1 , k 2 , k 3 ,pi,p 2 , a 1 ,a 2 
are determined by a photogrammetric procedure 
known as analytic self calibration 3 . This procedure 
is typically accomplished by bringing the recording 
camera into a calibration lab. The parameters 
w, <f, k, Xtj, Yq, Ztj are determined through a 
technique known as space resection. Measured 
control points, or targets, are imaged with the 
camera in its proper position for obtaining light 
sheet images. The resection technique involves 
solving equation 6 above for u>,<f>,K, Xq, Yq, Zq 
when both the image and object coordinates of the 
control points are known. 

The technique for determining the light sheet 
orientation parameters 6, differ for the experiments 
described below and no general analytic method has 
yet been incorporated. The light sheet parameters 
can be fixed by alignment procedures for simple 
geometries. In other cases transit measurements 
with respect to the control points must suffice. 

Until analytic techniques are incorporated, the light 
sheet parameters will limit the system accuracy. 

TRANSLATING LIGHT SHEET APPLICATION 

A wind tunnel laser light sheet flow visualiza- 
tion experiment was performed by researchers in- 
vestigating the relationship between buffeting and 
the vortical flow Held over a wing 4 . The investi- 
gation was conducted in the NASA Langley Re- 
search Center Basic Aerodynamic Research Tunnel 
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(BART). The BART is a test facility dedicated to 
code validation 5 . Off-body flow visualization was 
obtained from laser light sheets illuminating smoke 
entrained in the flow. The model, shown in figure 2, 
consisted of twin vertical tails placed on tail sup- 
port booms aft of a 76 ° flat delta wing. The model 
was tested at a = 20 ° and a = 30 ° . 

An illustration of the BART experimental setup 
for a = 20 ° is shown in figure 3. A stationary 
video camera mounted above the test section’s 
glass ceiling and forward of the model was used to 
record the experiment. The light sheet, mounted on 
a traversing mechanism above the test section, was 
directed downward onto the model. The beginning 
and ending light sheet locations are represented 
in the figure as grids. The 3D region over which 
flow Held data was obtained during the experiment 
is bounded by the camera field of view and the 
beginning and ending light sheet planes. The 
camera sensor array, represented by grid lines in 
the image plane, projects through the lens center to 
grid lines in the light sheet planes. 

In order to obtain the position and orientation 
of the recording camera via the space resection pro- 
cedure, a minimum of three control points are re- 
quired. These points were conveniently obtained 
for this experiment by using the three sharp cor- 
ners at the nose and wingtips of the model. The 
locations of these corners were obtained from model 
construction schematics and were assumed to be ac- 
curate. Figure 4a shows their locations in the video 
image. The normalized object coordinate system 
was defined with the origin at the model nose and 
the wingtips in the X — Y plane, as shown in fig- 
ure 5. The X— coordinates of the wingtips were 
defined to be 100 percent model chord with the Y — 
coordinates scaled appropriately. 

Light sheet plane location data is required for 
the reconstruction process. The task of accurate 
light sheet placement was simplified for this experi- 
ment by the flat triangular model shape. The sheet 
was manually aligned so that both the wingtips 
were in the light sheet and a flat mirror placed 
on the model surface would reflect the light sheet 
back to its source. This alignment procedure in- 
sured that the light sheet was parallel to the Y — Z 
plane. The traversing mechanism was used during 
the tests to position the light sheet along the A— 
axis from 20 to 160 percent model chord, in ten- 
percent increments. These positions are shown in 
figure 5. The equations for the light sheet planes 
are then given by A = 0.2, A = 0.3, . . . , A = 1.6. 
These equations correspond to equation 4 with 
&i = —1.0, &2 = 0,63 = 0, and 64 = 0.2, 0.3, . . . , 1.6. 

A typical light sheet image for this experiment is 


shown in figure 4b. While the geometric reconstruc- 
tion of the light sheet plane has been discussed in 
detail, the flow features revealed by scattered light 
and recorded as gray level variations in the image 
have not yet been addressed. The process of assign- 
ing a gray level value to each object point location 
based on the scattered light intensity of the corre- 
sponding image point will be described in the sec- 
tion on the computer visualization process. 

ROTATING LIGHT SHEET APPLICATION 

An in-flight flow visualization experiment was 
performed by researchers investigating the origina- 
tion of the vortex system and multiple vortices over 
the wing of the F-106B aircraft 6 . A photograph of 
the F-106B indicating the vapor screen components 
is shown in figure 6 . (The figure shows a remov- 
able leading-edge flap attached to the wing which 
was not used for the flight data presented in this 
paper.) Off-surface flow visualization was obtained 
by illuminating propylene-glycol vapor entrained in 
the flow with a rotating light sheet. The propylene- 
glycol vapor was introduced into the flow through 
a probe located underneath the port wing near the 
leading edge. The light sheet was emitted through 
a narrow rotating slit located on top of the fuse- 
lage. The rotating light sheet swept out a range of 
angles from 40 ° to 125 ° aft of the aircraft center- 
line, and remained perpendicular to the wing upper 
surface. The light sheet angle was recorded during 
sweeps. Two video cameras, one mounted on top 
of the fuselage in front of the light sheet housing, 
and the other mounted on the port inlet, were used 
to record the experiment. The inlet camera had a 
view of the apex region and the upper wing surface, 
while the top camera had a view more aft along the 
wing upper surface. 

An illustration of the F-106B experimental 
geometry for the top camera is shown in figure 7. 
This illustration is similar to the one shown for 
the BART experiment except that the camera is 
omitted for clarity. The 3D region over which flow 
field data was obtained is the intersection of the 
camera field of view and the area illuminated by 
the rotating light sheet. The object coordinate 
system shown was adopted to be compatible with 
other measurements performed on the aircraft. The 
techniques used to locate the control points and 
light sheet are described below. 

Since no conveniently identifiable wing features 
were available, such as the sharp corners on the 
BART delta wing, the control points had to be 
manufactured. A video image of the control points 
used is shown in figure 8 a. The control points were 
circular pieces of cardboard centered as close as 
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possible on the intersections of wing spars and ribs 
whose X and Y coordinates were documented in 
the original F-106B construction schematics. The Z 
values were measured by a scale from a level line. 
Given the luxury of a relaxed flight schedule, a 
photogrammetric survey of the aircraft would have 
been preferred. 

The light sheet was designed to rotate about 
a line parallel to the Z— axis. The measured lo- 
cation of the light sheet origin was (X s , Y s , Z s ) = 
(438.7,0.4, 51.5) in inches as shown in figure 9. The 
aircraft instrumentation system provided the rota- 
tion angle, <F, measured from the nose of the air- 
craft. The equation of the light sheet, as a function 
of $, is given by 

tan(q>) = ( 8 ) 
and equivalently in the form of Equation 4 as 
(tan(<S>))X + (-l)V + (0 )Z + (Y s - X s tan(<S > )) = 0. (9) 

The videotape recording of the in-flight ex- 
periment contains the primary data, i.e. , vapor- 
scattered light from the rotating light sheets. A 
typical image, representing one slice through the 
3 D flow Held, is shown in figure 8b. One factor that 
makes this image difficult to analyze is that the 
outline of the wing upper surface is barely visible. 
The light sheet footprint, i.e., the intersection of the 
light sheet with the wing upper surface, produced 
the most reflection and therefore appears as the 
brightest region in the images. The footprint pro- 
vides an indication of the light sheet position when 
examining the 2 D images. 

COMPUTER VISUALIZATION PROCESS 

The computer visualization process is a com- 
bination of digital image processing techniques for 
digitizing, enhancing, and analyzing sequences of 
wind tunnel and in-flight images recorded on video- 
tape, and interactive computer graphics for display- 
ing and analyzing the reconstructed image data. 

All visual data analysis was performed in the Data 
Visualization and Animation Lab (DVAL) at Lan- 
gley Research Center. The lab provides advanced 
visualization tools on an integrated system of high- 
performance graphics workstations and digital im- 
age processing equipment. DVAL resources include 
Silicon Graphics, Inc. (SGI) Iris workstations for 
performing interactive computer graphics, high- 
quality image rendering, and advanced animations; 
Sun workstations for performing image analysis; 
and a Video Image Processing System (VIPS) for 


performing video-rate image acquisition, processing, 
and storage. The VIPS is a combination of hard- 
ware and software designed to provide a special- 
ized capability for applying enhancements at video- 
rates to sequences of wind tunnel or in-flight images 
recorded on videotape. 

The overall computer process is illustrated in 
figure 10. Videotape frames were digitized and 
processed on the VIPS. The digital images were 
then analyzed on the Sun computer to extract 
vortex core locations. The images and the core 
locations were subsequently transferred to an Iris 
workstation where the triangulation software was 
executed to map image pixels into 3 D space. The 
camera and light sheet parameters are necessary 
inputs to the triangulation software. The Flow 
Analysis Software Toolkit (FAST) on the Iris was 
used to interactively display the 3 D image data 
along with the numerical surface geometry of the 
model. 

Digitizing and Processing the Videotape 

In the two examples described above, the vi- 
sual record of the light sheet flow visualization ex- 
periment is a standard S-VHS videotape. Prior to 
image reconstruction, these analog video signals 
must be converted to digital images. The VIPS, 
composed of an Imaging Technologies, Inc. (ITI) 
Series 150 modular image processing subsystem 
and a Storage Concepts, Inc. real-time digital disk 
(RTDD) subsystem interfaced to a Sun computer, 
was used to digitize, process, and store select frames 
from the raw video footage of the flow visualiza- 
tion experiments. The block diagram in figure 11 
illustrates the pipeline configuration of the modules 
within the Series 150 subsystem and the interface to 
the RTDD subsystem. The Analog-Digital Interface 
(ADI) converts the analog video signal into a two- 
dimensional array of numbers that can be passed 
sequentially through the pipeline to the other mod- 
ules in the subsystem. For each frame of video, the 
ADI divides 480 of the active horizontal video lines 
into 512 samples, known as pixels. Each pixel loca- 
tion is assigned an integer gray level value between 
0 and 255 corresponding to the brightness of the 
image at that location. The resulting 512 x 480 x 
8 bit digital images are passed through the pipeline 
to other modules in the pipeline where additional 
processing may be performed. The Real-Time Mod- 
ular Processor (RTMP) is the first module in the 
pipeline, followed by the Arithmetic Logic Unit 
(ALU), and finally the Frame Buffer (FB). Pixels 
are processed through the pipeline at video-rates, 
or 30 frames per second. The RTDD is capable of 
storing approximately eight minutes of full-frame 


5 



digital image data. The software interface is Noesis 
Vision, Inc’s VISILOG Image Processing Software. 

During the BART experiment, the light sheet 
was positioned at 15 different stations along the 
model, ranging from 20 percent model chord to 
160 percent model chord, as shown in figure 5. 
Approximately 20 seconds of video was recorded 
for each light sheet station. The ITI pipeline was 
configured to process the digital images coming 
from the ADI using a pointwise recursive digital 
filter of the form: 


9i = afi + (1 - (10) 

where /, and §i are the input and output images, 
respectively, and 0 < a < 1. This temporal digi- 
tal filter is commonly used in signal processing 7 to 
smooth digital data. A weight of a = 0.25 was se- 
lected with the goal of removing noise in the images 
without significantly altering temporal information. 
The first input image was /o and an image of all 
zeros was used as the initial value for g-\. The re- 
sulting time-filtered digital images were stored on 
the RTDD. The processed images were then played 
back from the RTDD at speeds less than video rates 
for inspection. The researcher selected one repre- 
sentative image for each light sheet station based 
on subjective image quality. The images shown in 
figure 4 are two examples of images digitized and 
processed on the ITI subsystem. 

For the F-106B in-flight experiment, the re- 
searcher previewed the videotape in advance, and 
selected specific frames to be digitized. Frames of 
interest were identified by the visual time code ap- 
pearing on each video frame. No processing was 
performed on these images prior to image recon- 
struction. 

Extracting Vortex Core Locations in 2 D Images 


positions (m 2 , n 2 ), . . . , {mjg, «jv) by 
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( 11 ) 


Since the center of a dark region is of interest, the 
gray level values of the pixels in this region are 
contrast reversed, i.e. , replaced with 255 — for the 
centroid computation. The m ^ and n c f, coordinates 
are the vortex core centroid. 

The center of brightness computation is re- 
stricted to a subset of pixels comprising the dark 
interior region of the vortex. This subset is deter- 
mined using a sequence of simple steps. First, the 
researcher specifies an initial estimate of the vortex 
core location by selecting a single pixel somewhere 
inside the dark interior region of the vortex. This 
initial guess pixel serves as the input to a technique 
which finds the valley of the region surrounding this 
pixel, analogous to letting a ball roll down a hill. 
The output pixel location replaces the researcher’s 
initial guess and serves as the input to a technique 
which searchs for the ridge of the hill. The result is 
a closed boundary, shown in figure 12b, defining the 
subset, A, of pixels to be considered for use in the 
center of brightness computation. Pseudocolor, an 
assignment of colors to gray level values, is used in 
the figure to emphasize contrast variations within 
the vortex core. Once A has been determined, B, 
a subset of A, is defined by applying a gray level 
threshold to the pixels in A. Subset B, shown in 
Figure 12c as an entirely black region, represents 
the pixels in subset A with a gray level value below 
the chosen threshold. The pixels in B are used to 
compute m c 5 and n c f, from the equations above, af- 
ter contrast reversal. The crosshair in Figure 12c 
indicates the computed pixel location of the vortex 
core. 


Once the light sheet images have been digitized, 
interactive image analysis techniques can be applied 
in order to extract quantitative information such 
as vortex core locations. The BART light sheet 
image at station X = 1.1 will be used as an 
example. An enlarged gray scale image of the port 
side vortex is shown in figure 12a. The vortex core 
appears as a dark interior region, void of seeding 
particles. The location of the vortex core centroid 
is determined by finding the center of mass of 
the pixels within these dark regions. A better 
terminology might be center of brightness, since 
a brightness value is more commonly associated 
with pixels. The horizontal scan line and column 
position of the center of brightness is denoted by 
(m c 5 , n c j). This coordinate is determined from N 
pixels having gray level values /i , I 2 , ■ ■ . , In, and 


Interactive Display of 3 D Data 

The photogrammetric reconstruction of one 
512 x 480 digital image yields a 3 D coordinate 
and a gray level value for every pixel, requiring 7.5 
Mbytes of storage space. The total storage require- 
ment depends on the number of light sheet planes 
to be analyzed and the additional storage required 
for the numerical surface geometry of the model. 
Typically, anywhere between 10 and 15 images per 
test condition were studied for the two experiments 
discussed above. A sophisticated visualization pack- 
age is required to interactively display experimental 
data sets of this magnitude. The package chosen 
for this work is the Flow Analysis Software Toolkit 
(FAST) running on the Iris workstation. FAST was 
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originally developed as a post-processing graphics 
environment for viewing and analyzing large 3 D 
computational fluid dynamics (CFD) databases 8 . 
The package allows researchers to load in CFD 
grids and solution data sets stored on the Iris in 
PLOT3D format 9 , and to visualize this data. In 
this work, experimental data is stored much the 
same way CFD data sets are. A data set containing 
the 3 D coordinates of object points corresponding 
to video image points is output by the triangula- 
tion software. This data set is also a PLOT3D grid, 
and can therefore be loaded into FAST in the same 
manner a CFD grid is loaded. A data set contain- 
ing image intensities is output during the digiti- 
zation process. This data set can be thought of as 
containing the functions values for each grid loca- 
tion, and can also be loaded into FAST in much the 
same way a CFD solution data set is loaded. 

There are many features of the highly interac- 
tive FAST environment that are useful for analyzing 
the experimental data sets. Figure 13 shows sev- 
eral user interface windows that appear on the Iris 
display during a FAST session. The delta wing nu- 
merical surface geometry, reconstructed light sheet 
images, and vortex core paths are shown in the Ac- 
tive View Window. In order to create this display, 
grid Hies for the delta wing, twin tails, light sheet 
stations, and vortex core locations are selected for 
data input. Also, the light sheet intensities corre- 
sponding to the light sheet grid are selected. Pseu- 
docolor is used to display the light sheet intensi- 
ties. A range of low intensities, (from 0 to 20) are 
clipped, i.e. , not displayed, using a thresholding op- 
tion in FAST. This allows only the vortex struc- 
ture to be seen. Also, object points with Z values 
less than 0 are not displayed, eliminating unimpor- 
tant data below the wing. The entire scene can be 
rotated, translated, and zoomed in and out using 
the interactive viewing controls. A data probe and 
probe typeout can be used to examine the 3 D coor- 
dinates of any grid point in any of the scene objects. 
In the figure, the probe is being used to extract the 
X, Y, Z value for the port side vortex core at light 
sheet station X = 0.7. 

The entire scene shown in the Active View 
Window in figure 13 is created in a progression, as 
shown in figure 14a-f. First, the numerical surface 
geometry of the delta wing and left and right tails 
is read in and displayed with grid lines (figure 14a). 
The model can also be rendered as a smooth shaded 
surface to provide a more realistic view (figure 
14b). Figure 14c shows the smooth shaded surface 
with the light sheet grid for station X = 1.1. 

The function values at each grid point, in this 
case image intensities, are interpolated to produce 


a smooth shaded surface as shown in figure 14d. 
Function values below a specified value are clipped 
and every other light sheet station is displayed in 
figure 14e. Lastly, figure 14f shows a continuous 
vortex core path by connecting the vortex core 
points extracted for each light sheet station with 
straight lines. A similar progression for the F- 
106B port wing is shown in figure 15a-f. The entire 
F-106B aircraft can be displayed in FAST with 
selected light sheet images as shown in figure 16. 
Details of the F-106B aircraft, such as the nose, the 
fuselage, the light sheet housing, and the light sheet 
are accentuated in this view. 

The objects in a scene can be rotated and 
viewed from any angle. Figure 17 shows a top and 
side view of the delta wing with just the recon- 
structed vortex core points, connected by straight 
lines. The top view gives insight into the distance 
the core path is from the leading edge, whereas the 
side view provides insight into the height of the vor- 
tex core above the delta wing surface. 

An advantage of using FAST to analyze the 
experimental data sets is the capability to perform 
visual comparisons with CFD data sets. Figure 
18 shows the delta wing with experimental data 
on the port side and computational stagnation 
pressure contours on the starboard side 10 . The 
computed yellow particle traces indicating the path 
of the vortex core show good correlation with the 
experimental results. 

CONCLUSIONS 

It has been demonstrated that the computer- 
aided image reconstruction technique can be effec- 
tively applied to research projects using moving 
light sheets to visualize 3 D flow fields. The tech- 
nique provides researchers with a new 3 D quanti- 
tative analysis capability for 2D light sheet images 
recorded on videotape that were previously diffi- 
cult to evaluate. Once represented in 3 D, the im- 
ages can be observed from any viewpoint, freeing 
the researcher from the original camera view. In 
addition to a powerful 3 D visualization tool, the 
technique provides mathematically precise quanti- 
tative information about locations and structures 
of vortex systems. FAST allows the X, Y, Z loca- 
tions and function values of image points and model 
grid points to be examined interactively. Informa- 
tion, such as the distance from vortex core locations 
to the wing leading edge, is easily extracted. This 
technique also shows promise in aiding the visual 
comparison of experimental and CFD data. Exper- 
imental data is visualized using the same computer 
graphics environment, FAST, commonly used for 
analyzing and interpreting CFD data. 
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